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Starting from the quasi one-dimensional kinetic energy 
of the d yz and d zx bands we derive a bosonized description 
of the correlated electron system in S^RuCU. At interme- 
diate coupling the magnetic correlations have a quasi one- 
dimensional component along the diagonals of the basal plane 
of the tetragonal unit cell that accounts for the observed neu- 
tron scattering results. Together with two-dimensional corre- 
lations the model consistently accounts for the normal phase 
specific heat, cyclotron mass enhancement, static susceptibil- 
ity, and Wilson ratio and implies an anomalous high temper- 
ature resistivity. 



I. INTRODUCTION 

Sr2Ru04 is the first layered transition metal oxide thai, 
exhibits superconductivity in the absence of copper ions.EJ 
The lattice symmetry is tetragonal and isostructural to 
La2CuC>4 with lattice parameters a = b = 3.87 A in the 
RuC>2 plane and c = 12.74 A out-of-plane. No structural 
instabilities are-pbserved.B The first de Haas-van Alphcn 
(dHvA) resultsd and band structure calculations in local 
density approximation (LDA)o show three bands cutting 
the Fermi level with quasi two-dimensional Fermi sur- 
faces. They can be mainhi-associated with the three ti g 
orbitals of the Ru 4+ ions,ErB and are consistent with the 
metallic properties and the strongly anisotropic transport 
along the c axis .El 

The enhanced specific heat and magnetic susceptibility 
indicate the presence of significant correlations .El Consis- 
tently, results from, angle resolved photo emission spec- 
troscopy (ARPES)lJ and dHvA measurementsa suggest a 
strong electronic mass renormalization. The material is 
FemiLliquid like in a temperature range of T c < T < 30 

The significant correlations in Sr2Ru04_.the 5=1 
moments on Ru 4+ impurities in S^IrO^,!! 1 ! and ferrp* 
magnetic correlations in SrRuC>3 led Rice and SigristE3 
to propose that the superconducting order parameter 



has p-wave symmetry promoted by ferromagnetic corre- 
lations analogous to 3 He. The absence of a-daange in the 
Knight shift in the superconducting phaseEjO supported 
that notion as well .as the temperature independent mag- 
netic susceptibilitytJ and the enhanced relaxation time 
in Muon spin resonance (^SR)N at T < T c ~ 1.5 K. A 
similar proposal was made by Baskaran based on a com- 
parison with high T c materials and emphasizing the role 
of Hund's rule coupling.lill 

Since then tremendous experimental effort has been 
made trying to verify the predicted p-wave symme- 
try of tly ..superconducting order parameter. Neither 
tunnclingE3o nor thermal conductivity experimentscjo 
or ac-susccptibility measurements^ under different mag- 
netic field geometries gave conclusive proof of the anal- 
ogy to 3 He. No indication fopJew-omaguelie correlations 
has beenJauiid in ARPESjHSEa LDA,£§El or neutron 
scattaiMjSE] investigations. Furthermore, the sgwscific 
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heatE2rE3 and nuclear quadrupole resonance arc 
consistent with two-dimensional gapless fluctuations in 
the superconducting phase of Si'2Ru04 which are absent 
in superfluid 3 He. 

The controversy about the proper description of the 
electronic correlations in Sr2Ru04 is reflected most im- 
pressively by the variation of values of the on-site 
Coulomb repulsion U used in the mostly perturbative 
approaches to match experimental results. Examples 
U « 0.42-eV^3 n_2 eV,H 1J2^ 1.5 eV,B 0.345 eVP 
0.175 eV,E3 2 evj?il 0.048 eV.H Comparing these values 
of the interaction with the bare Fermi velocity of vf « 0_Z 
eVa from band structure calculationscH and ARPESE3 
points toward an intermediate coupling regime. 

In contrast to the effects of the interactions the 
bare electronic band structure, has beep-, determined 
unambiguously from dHvAjilrej ARPES,cfl and x-ray- 
absorption pieasurementsB in consistency with LDA 
calculations .oE3 The overlap of the electronic wave func- 
tions of the—rL^-and d yz orbitals is dominantly one 
dimensional Ej'lZIQ The interaction and additional hop- 
ping channels lead the d zx and d yz electrons to hybridized 
into two bands. Their Fermi surfaces are referred to as 
the a and sheets. The electwwj* in the d xy orbital form 
the two-dimensional 7 sheet .13123 

Correlations in effective one-dimensional systems show 
power law behavior E3 They are always more singular 
than two-dimensional correlations which diverge at most 
logarithmically.! 4 -!! Since the kinetic energy of the d zx and 
d yz electrons is quasi one-dimensional we expect their 
correlations to play a dominant role. 
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The quasi one-dimensional kinetic energy of the 
d zx and d yz electrons allows for the bosonized, non- 
perturbative description of the low energy electronic ex- 
citations. This description is introduced in Sec. |ll| a nd 
its fundamental properties are discussed. Section is 
devoted to the expected corrections from hybridization 
effects and the 7 sheet that have been neglected in the 
initial model. The comparison with experimental results 



in Sec. IV reveals the qualitative and quantitative consis- 



tency of the model within the framework of its applicabil- 
ity. A comparative discussio n of a lternative perturbative 
approaches is included [Sec. IV B]. 

The present paper is part I of a series of three. Part II 
(Ref. |^|) is devoted to the superconducting phase. The 
in-plane correlations are described via the model derived 
herein. The inter-plane pair-correlations are enhanced as 
a consequence of the body centered crystal structure and 
can be treated mean-field like. 

Part III (Ref. ^) consistently explains the experimen- 
tally observed unconventional transitions under magnetic 
fields based on the model derived here and in part II. 



II. SUBSYSTEM OF d zx AND d yz BANDS 

The band structure as determined from .dHvA,0 
ARPES,E-3 and x-ray^ibsorption measurementsO as well 
as LDA calculationsE3'E3 together with the anticipated 
intermediate interactions suggest a three band Hubbard 
Hamiltonian as the generic model. 
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Here cr^* z ' denote the 
^ 2 and 2L 



Pauli matrices with cr° a , — 
(°a a') 2 an d 2L = L x = L y is the linear dimension 
of the system. We limit the description here to the 
RuC>2 planes and generalize when necessary to a three- 
dimensional array of planes. The intra-orbital Coulomb 
repulsion is larger than the inter-orbital repulsion, i.e., 



Jjv=v 



U 



Ulit and U > U 2 



U, 



> Ul-u^, 

Hund's rule coupling lowers the inter-orbital Coulomb 
repulsion for electrons in a spin-triplet configuration with 
respect to the spin singlets. The full treatment of the 
involved exchange interaction terms within the frame- 
work of the bosonization approach discussed in Sec. [I A 



is rather involved and only possible in approximations. A 
qualitative study of the effect of Hund's rule coupling is 
possible by setting U\ > U 2 and neglecting the exchange 
terms. The expected corrections due to ex change and 
other terms are discussed closer in Sec. |II C| . 

The Hamiltonian Eq. (g) is SU(2) invariant both in the 
spin sector and in the orbital sector yielding an effective 
SU(2)(g)SU(2) symmetry. The orbital degrees-etireedom 
are sometimes referred to as electron flavors. tiZrliiJ 



A. Bosonization 



In this notation the electron creation and annihilation 
operators are cj v and q v for orbital v with spin a on 

site Z, ni M ^ a is the usual electronic density operator, t^, 
is the hopping matrix element between site Z and Z', and 
U v ' V i is the on-site Coulomb repulsion. 

As discussed in the introduction we expect the inter- 
esting low temperature physics to be dominated by the 
quasi one-dimensional d zx [y = x) and d yz (y = y) 
bands. We retain here only the dominant hopping am- 
plitudes to = ti'i + ~ 

hybridization of the bands later in Sec 



tf'i+y an d discuss effects from the 
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The con- 
tinuum representation is introduced via q,„ i(T — * tp u , a (r) 
with p u ,a{f) — 4>t a { r )i } v a-( r )- The bands are linearized 
with Fermi velocity v-p w \/3 to- Note that the "veloci- 
ties" in the present paper define the kinetic energy scales, 
i.e., v-p = vp/a (h= 1). 



In Eq. (g) the d zx fields only have a kinetic energy 
contribution along the x direction while the d yz fields 
only have a kinetic energy .along y. We can thus switch 
to the chiral representationc3 



W,.(r) = iW)e lfc ^ +V(r), 



(3) 



for each orbital degree of freedom v at each position 
r v '^ v transverse to the propagation. The right (R v . a ) 
and left (L^ ^.) moving fermions of each species can now 
be bosonized.cZrCJ 



Ru,a{r) = lim 



^A Tv ') „-WH[6 v , a {r)+4> v ,a{r)\ 



(4) 



L u , a (r) = lim e -i>/S[«,.,M-^W] (5) 



a— »0 

Here (f>v,rr( r ) are the Bose fields with their conjugate mo- 
menta H v<(7 (r) = d v Q v . a {r) which satisfy the commuta- 
tion relation 

[<f>»Ar),~n-»>,<r'(r')]=i6 v y5 a ^8(r-r'). (6) 

The Klein factorsE3 ?j„,a-(?v) assure the proper commu- 
tation relation between the different fermion species, a is 
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a short range cut off associated with the in-plane lattice 
constant. The bosonized Hamiltonian can be written as 



H *> = I™ E E AM? (n^ + ^-M 2 

t^oo ct^ct' „^„' •'-i L 

, U ° l\a a ir« ^ 1 , cos \/47r(0«/,<7 - </v,i7') 
+ — I [d v (j> v ^\[d v (j) v ^\ + 



(27ra) 2 



■[dv^>v,(t][dv'<t>v> 



4(27ra) 2 

: f cos [V47r(0„ iO . - i<r / ) - 2/c F (r„ - r^/ )] 
+ cos [V47r(^, CT + ^r/', CT ') - 2k F (r u + rv)]J 



+ -Y[du4>u,o][dw(t>u',o] + A( ^ 2 , 2 
4 4(27ra)' ! 



(cOS [V^TT{4>y.a ~ 4>v\<j) ~ 2ks(r v ~ TV')] 

+ cos [a/47t(^ (T + 0„' )O -) - 2fc F (r„ + rv)] 



(7) 



The standard approach to separate spin and charge 
degrees of freedom in Eq. (Q) is to introduce charge and 
magnetic fields for each flavor via 



fu,m{r) = [<i>v,i(r) - ^(r)]/V2, 
<Pv,c(r) = [^,T(^) + ^,l( r )]/V2, 



(8) 
(9) 



respectively. TL Utm and n„ c are the corresponding con- 
jugate momenta. The bilinear part of the Hamiltonian 
H2D is composed of the charge and (magnetic) part 



H, 



1 



c(m) 



E 



d 2 r 



v F Il 2 



y,c(m) 



E ) v }y ( 9 »' f>w ,0(11: i 



) ■ (10) 



The matrix elements for the charge (magnetic) part are 

= VF + (_)[/•„ and V$ m)) = 
Ui + (— )Ui, Equation (10) is the Hamiltonian 



given by Vx% 

<0) 



of a crossed sliding Luttinger liquid studied in Refs. |5l| 
and ^2| The authors find no significant change in the 
decay of the low energy correlations with respect to the 



one-dimensional case where V a 



(c(m)) 



■r,y 



0. A perturba- 
tive treatment suggests that the inclusion of the inter- 
action term in Eq. (Q) leads to two-dimensional correla- 
tions which still decay algebraicallycil In the absence of 
Hund's rule coupling, where U2 — U±, the magnetic sec- 
tor fields ip Xjm and ip Vt m are decoupled. Similar models 
are obtained for coupled Luttinger liquids .EjEj 

To study the qualitative properties of the model de- 
fined by Eq. (Q) with parameters relevant for Sr2Ru04 it 
proves useful to use the symmetry of the orbital degrees 
of freedom. We introduce charge (/1 = p), spin (/i = s), 
flavor (/1 = f), and spin- flavor (/z = sf) fields via the 
canonical transformation 



M r ) = \ E a l,v° h o.a <t>»A r ) > 



(11) 
(12) 



The matrices are (a, 6) = {(0, 0); (0, z); (z, 0); (z, z)} for 
(1 = {p;s;f;sf}, respectively. The fields M (r) are identi- 
cal to those of the resonant-levelrjaaodel used to describe 
the two-channel Kondo problem.cH The fields are simple 
linear combinations of the charge and magnetic fields, 
e.g., 



(j) s (r) = (ip Xtia + <p y ,m) /V2 
(f> s {(r) = (<p x>m - <Py,m) /V2. 



(13) 
(14) 



Note also that the charge and spin sector fields are sym- 
metry related via the reflection y — > — y R yM <-> 
L y ,a 4>i <-* <i> P and <p s <-> <j> s f. 

The representation can be simplified by introducing 
the variables x — (x + y) and y = -^(x — y) with 

r = (x, y~y. The charge Hamiltonian in Eq. ( JlCj ) becomes 



Hc = \ J d2r { VF ( H P + n 



V c [d x ^> p + d v ct> { } 2 + 1/ c [%0 P + cV0f] 2 }, (15) 



while the magnetic Hamiltonian is 
1 



-Id 2 , 



{v F (n 2 s + K 



+ V m [dxfis + dy(j) si ] 2 + V m [dy4> s + cV^sf] 2 } ■ (16) 



The energies are 



V c = v F + U + {U 1 + U 2 ), 
V c = v F + Uo-(U 1 + U 2 ), 
V m = v F - U + {Ui - U 2 ), 
V ra = v F -U -{U l -U 2 ). 



(17) 
(18) 
(19) 
(20) 



Applying Eqs. (pT) and (|lj) the interaction term in the 
Hamiltonian Eq. (% factorizes into contributions of the 
four spin, charge, flavor, and spin-flavor degrees of free- 
dom. 



H 



U 



int 



d r cos v A-K<j) s {r) cos v 47r</> s f (r ) 



(27ra) 2 

j- 1 ,„ / d 2 r cos \V^TT(f>{ (r) - 2V2k F y] 
(2nay J 

x(Ui cos V4it(j) s (r) + U 2 cos %/4~7n/> s f (r)^j 
j^- 2 J d 2 r cos [V4^> p (r) - 2V2fc F ^] 
x (1J2 cos V^(t)s(r) + U\ cos \/~An<i> s i (r)^ 



(21) 
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The limit a — > and L 



Hamiltonian of the d zx -d yz 



density wave instabilities are ob- 
- s Tlic values for t ne,og-g^^uloiii! i 



oo is understood. The total 
subsystem is i?2D — H c + 

H m + -f/jnt • 

No spin or chars 
served in Si^RuCi. 

repulsions discussed in the literaturdiJ123'E3'ti3£!i] point 
toward an intermediate coupling regime if compared to 
the bare Fermi velocity of vp « 0.7 eV from band struc- 
ture calculations.cH For repulsive interactions the oper- 
ators in Eq. have been shown to be-jnarginally ir- 
relevant both in one and two dimensions J13 The physical 
properties are therefore determined by H c and H m with 
quantitative corrections from flint- Corrections from hy- 
bridization terms not included in the bosonized model 
are discussed in Sec. 



The Luttinger liquid parameter K e g and-the velocity 
v e g are effective parameters of the theoryEll but can be 
associated with 



K. 



and 



cir 



Vctt 



VVm/VF 



V F . 



(24) 



(25) 



III 



Note that Eqs 
invariant under the reflection x *— 
to x -> x, y -> -y, 



15|) through ( pl| ) are still explicitly 
y which is equivalent 
and 4> s f — ► — (j) s {. The 



9f 



same applies for 
4>f <-> 4> n and <b R * 



-y where x 



y- 



FL, 



Li, 



B. Effective one-dimensional model 

In the intermediate coupling regime the model defined 
by Eqs. (|l5|) through (|2l|) e xhibits a number of singular 
points for V c ( m ) = or V c / m ) = 0. For Uq < vp and 
Uq > U\ > U2 the relevant limit is V m — > 0. Then the 
magnetic Hamiltonian Eq. (|l6|) only has terms in c^-</> s 
and dy<p s {- 



1 



d 2 r 



{vf (n s 2 



The model Eq. (|22j) is explicitly invariant under the 
transformation y — > —y, where x «-> y and (f> s <-> <f> s f. The 
effective one-dimensional model Eq. (^3|) describes the 
quasi one-dimensional magnetic correlations along both 
in-plane diagonals of the tetragonal unit cell. 

In realistic systems V m > and the effective model is 
applicable only for sufficiently large temperatures T > 
V m . A close discussion of the values appropriate for 
Sr 2 Ru04 is given in Sec. IV. 

The total Hamiltonian of the low energy in-plane d zx - 
d yz correlations is given by H2D = H c + H c g + iJ; nt in 
Eqs. (|l5|), (p3[), and (|2~l|). The interaction term H int can 
be neglected when appropriately rescaling the parameters 
K cS - K* cS , V c -» V*, and V c -» V* c . 

After eliminating the term ~ {dx4> s ){dy4>sf) in Eq. ( p^ ) 
by a sliding transformation the magnetic Hamiltonian 
can be written as a superposition H m » H R + H s { w H e g 
where H s and H s { are obtained from Eq. ( p3|) by replacing 
(eff — > s) and (eff — > sf , a; — > y), respectively. Note that 
here v s = u s f = u e fif/2. In the case of SU(2) symmetry in 
the spin subspace the interaction term Eq. ( pl[ ) yields a 

In Rcf. 



rescaled K,, 



n s 2 f ) 



K* = 



1. 



45 it is suggested that 



the spins of the electrons in the d zx -d yz subsystem are 
in an easy-plane configuration which implies K* < 1 and 



(22) 



The representa tion of the spin and spin-flavor fields in- 
troduced in Sec. II A has the property that along x = y 
one finds dy4> s = dx4>si = 0. This becomes obviousJrom 
Eqs. (13) and (|l4| ) together with the symmetrycj im- 
plied relation d y ip y m \ = d x cp x m \ . Thus, the fields 
4> s and (f) s { indeed depend only on y and x, respectively, 
as implied by Eq. (|22| ) and H m is effectively one dimen- 
sional. Since the spin and the spin-flavor channel are 
symmetry related the one-dimensional correlations in Eq. 
(E2n can be effectively described by 



ffcff = 



VctfL 



J dx[K c& Iil s + K^{^ cS f 



(23) 



Note that the coupling term {dx4> s ){dy<j>si) in Eq- ( HI ) 
can be eliminated by a mean-field decoupling and sub- 
sequent sliding transformation!^ This approximation is 
less severe than it may appear at first sight since along 
x = y the x dependence of </> s f and the y dependence of 
4> s can be neglected. Consequently, along x — y, (dx4>s) 
is just a y-independent constant with respect to (dy(f> s {) 
and vice versa. 



k; { < 1. 



C. Comment on Hund's rule coupling 

The treatment of the full SU(2) invariant Hund's rule 
coupling termEZI JhS^Sv'^i/ is difficult in the bosonized 
model. The presence of SU(2) symmetry breaking 
Dzyaloshinskii-Moriya interactions is obvious npffi the 
observed anisotropies of the static susceptibilitiesEjEI and 
is consistent with the expected presence of spin-orbit cou- 
pling in the Ru Ad orbitals.E3Hc3 Additional corrections 
are expected frcpa Kaplan Shckhtman-Entin-Wohlman- 
Aharony terms.EJ 

The crucial physical implication of Hund's rule cou- 
pling is that it couples the magnetic degrees of freedo m of 
the different orbitals. The model introduced in Sec. HA 



incorporates this effect qualitatively as becomes appar- 
ent from Eq. @. Notably is V m ^ V m only for Ui ^ U 2 
and consequently can the quasi one-dimensional model 
discussed in Sec. II B only be found in the presence of 
Hund's rule coupling. The quantitative justification of 



the model is shown phcnomcnologically in Sec. IV and in 
Refs. EBI and EI 
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Note rttat the estimated value for Hund's rule 
couplingEl in Sr 2 Ru0 4 of J H « 0.2 - 0.4 is larger 
than the estimate for the spin-orbit couplingpl of A « 0.1 
eV. Consequently the effective model for the magnetic 
correlations derived here is applicable even in the pres- 
ence corrections from spin-orbit coupling that lifts the 
degeneracy of the d zx and d yz orbitals because the larger 
Hund's rule coupling overcompensates the effect. 



III. INTER-PLANE COUPLING AND BAND 
HYBRIDIZATION 

The model described in Sec. |ll| has been based on 
the two d zx and d yz bands that are coupled via the on- 
site interaction only. In this section we discuss the d xy 
band as well as the d zx -d yz inter-chain and hybridization 
terms which are expected to qualitatively change the low 
temperature physics of the bosonized model. Here we 
discuss the magnitude of the terms and their expected 
impact. In Sec. |y|wc then discuss at what temperature 
which properties of Sr2Ru04 are determined by a certain 
subsystem. 

The in-plane resistivity is two to three orders of magni- 
tude smaller than that along the c axis.cl Consistently the 
dispersion of the Fermi energy along c is about 1% of tha 
in plane dispersion as probed be dHvA measurements. til 
Band structure calculations lead to an estimated hit 
plane hopping of about 10% of the in-plane hopping.l 

The appropriate inter-plane Hamiltonian with hopping 
amplitude t± is 



£ C l,v,<J C V ,v,a ' 



(26) 



V,v'= X , y l,V ,<T 



with only nearest neighbors Ri> — Ri + ^(±a, ±a, icy. 
The inter-plane hopping of the d xy band is an order of 
magnitude smallercil as a gonsequence of the in-plane ge- 
ometry of the d xy orbitalsB and can be neglected. Fourier 
transforming the Fermi operators via 



1 



leads to 
H± = 8t±_ £ 1 



(27) 



(28) 



The in-plane kinetic energies of the d zx and d yz electrons 
are 



(29) 



The total Hamiltonian for the d zx and d yz electrons 
H x + H y + H± can readily be diagonalized and yields 
the dispersion for the a and (3 bands as 



= E + t (cos ak x + cos ak y ) + 8t±gi_k 



The abbreviations <?o,fc = cos(afca;) — cos(ak y ) and g\ t k = 

cos 2 | a cos 2^ cos ^f- were introduced for lucidity. 

Panel (a) of Fig. |l| shows the resulting tight binding 
bands for k z — as a function of k x and k y . The param- 
eters are E = 0.22 eV, t = -0.3 eV, and t±_ = -0.02 
eV. The dispersion of the d xy or 7 band is given by 



E, 



(7) 



eV 



-0.39 — 0.54(cos ak x + cos ak y ) 



-0.44 cos ak x cos ak y . (31) 



cos ak x cos ak y stems from in-plane next 



The term 

nearest neighbor hopping. The corresponding real space 
hopping parameters are i 7 j 7 — —0.39 eV, i 7 i+s = 
tf'iir, = -0-27 eV, and tl'?!, , - = -0.11 eV. The qual- 
itative agreement with ARPES,Ea LDAP and dHvAQtl 
results is satisfactory. 

In order to obtain a more precise match with the 
three-dimensional Fermi surface suggested by dHvA 
measurementsa'cll it is necessary to extend the dispersion 
for the a and f3 sheets in Eq. (|3(]). The non- vanishing 
coefficient fco2 in Ref. |4^ suggest a term 2t z cos ck z 
from next nearest layer hopping. Spin-orbit coupling in 
the d zx - d^ .subsystem leads to a hybridization of the 
orbitals lSEZ J The fields in the Hamiltonians discussed 
in Sec. II A are describing hybridized d zx and gL_ or- 
bitals as a consequence of the on-site inter action. Lj In 
the framework of the tight binding model these effects 
can be modeled by introducing the on-site hybridiza- 
tion th- The in-plane dispersion is also enhanced and 
can be modeled by extending the diagonal contribu- 
tions Eq. ( p9| ) to include an effective inter-chain hopping 
to cos ak v — > to cos ak v + ti cos ak v i^ v . The resulting dis- 
persions are 

^ = + {to + ti) (cos ak x + cos ak y ) + 2t z cos ck z 
1 



(32) 



8t x g hk ± - V^5o, fe + 256ii ff f ifc + Uj L 



An appropriate choice of parameters is Eq = —0.29 eV, 
to = 0.3 eV, U = 0.03 eV, t z = 0.02 eV, t± = 0.02 eV, 
and th = 0.06 eV. The resulting bands are shown in Fig. 
panel (b) for k z = 0. 

The tight binding analysis leads to the following con- 
clusions. 

(i) The two- and three-dimensional corrections to the 
quasi one-dimensional d zx and d yz bands are of the order 
of 10% or 0.03 eV leading to the presence of the a and 
f3 sheets of the Fermi surface. Luttinger liquid behavior 
should only be observable at sufficiently high tempera- 
tures. For T > ~ 400 K the out-of-plane trans- 
port is incoherent and saturates while the in-plane resis- 
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(a) q?\ '<■ jr'j 3 




^5 



(b) /' 5 



r?5* J -0.6 !S 




FIG. 1. Tight binding model for the three bands that form 
the Fermi surface, (a) Reduced model with d zx -d yz hybridiza- 
tion trough inter-plane coupling Eq. J^ ). (b) Model including 
next nearest layer hopping, inter-chain coupling, and interac- 



tion induce on-site d zx -d yz hybridization Eq. 



for k z = 0. 



tivity is determined by the crossed sliding Luttinger liq- 
uid with linear temperature dependence. Consequently 
the model is consistent with the observed anomalous 
high temperature resistivityEa Also consistent is that the 
quasi-particterpeaks observed in ARPES disappear above 
T ~ 160 kHEI 

(ii) Since the on-site interactions are an order of mag- 
nitude larger than the tight binding parameters, i.e., 
Uq i.2 S> ti h _l z , dominant correlation effects are still de- 



termined or at least influenced by the Hamiltonians (23), 
(|l5|), and ( pl| ) with properly renormalized parameters. 
Examples are the mag netic s tructure factor (Sec. IV A ), 
the specific heat (Sec. IV C 2) and the degenerate super- 
conducting saddle point discussed in Ref. |45j. An account 
of the temperature dependence of the corrections due to 
the hybridization terms ti t h,±,z is given at the beginning 
of Sec. 0. 

(iii) As discussed in Ref. ^ the inter-plane hopping 
is the important parameter for the mean field supercon- 
ducting transition and is estimated to be t± ~ 20 meV. 
The resulting coefficient tj_/vp ~ 6 K of the inter-plane 
pair hopping term is consistent with the transition tem- 



perature of T c = 1.5 K on the mean field level. 



IV. APPLICATION TO Sr 2 Ru0 4 

The model derived in Sec. || readily accounts for the 
observed normal phase properties in S^RuG^. The most 
striking qualitative evidence is the recently discojtered 
scale invariance of the magnetic structure factored in 
agreement with the implications from the effective one- 
dimensional model Eq. (p^). The scale invariance has 
been observed for T = 60, 110, 160 K while it starts to 
break downed at-ff-j= 10 K which is in the regime of Fermi 
liquid behaviorJaeJ 

Therefore the hierarchy of the applicability of the 
model derived in Sec. |n| can be summarized as follows. 
For T > 400JC-we expect crossed sliding Luttinger liq- 
uid behavior .EiE3 Curvature corrections to the linearized 
bands of the order of 10% are conceivable at T ~ 700 
K. For 400 K > T > 25 K the system gradually crosses 
over to the Fermi liquid regime because the various cou- 
pling terms discussed in Sec. Ill become relevant at differ- 
ent temperatures.E-3 The observed scale invariance of the 
magnetic excitations suggest that these relevant terms 
mostly impact the quasi-two dimensional charge channel 
given by Eq. |l5|. The one- to two-dimensional crossover 
of the magnetic subsystem given by Eq. |l6| is determined 
by V m . Since the charge and magnetic channels are 
coupled via Eq. ( ]2~l] ) the Fermi liquid behavior only is 
fully observed in the electronic channel when the quasi 
one-dimensional magnetic fluctuations are frozen out for 
T V m . From the^£xperimentally observed onset of 
Fermi liquid behavioialid we estimate V m sa 25 K. Equiv- 
alently, a crossover to non-Fermi liquid behavior on en- 
ergy scales w > 2 meV is expected. 



A. Incommensurate magnetism 

An important probe for the interaction effects in cor- 
related electron systems is the magnetic structure fac- 
tor determined by neutron scattering. The quasi one- 
dimensional model derived in Sec. II B accounts for the 



dominant features of the magnetic response. 

The bosonization approach-correctly describes excita- 
tions near the Fermi surface pa i.e., for momentum trans- 
fer q ~ and q ~ 2fcp. The relevant momentum transfer 
for antiferromagnetic magnetic excitations^ is q ~ 2fcp 
since the integrated intensity of the structure factor for 
q — > vanishes E3 A one-dimensional model analogous 
to that of Eq. (H|) with equivalent bosonized Hamilto- 
nian and incommensurate back scattering wave vector 
q — is t' ie antiferromagnetic Heisenberg chain in a uni- 
form fieldO We thus expect the effective one-dimensional 
Hamiltonian Eq. ( p2| ) to describe an excitation spectrum 
as determined in Ref. [38] and sketched in Fig. ^. Since 
Eq. (^||) is a one-dimensional model along the diagonal 
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FIG. 2. Sketch of the spectrum of the elementary mag- 
netic excitations of a Heisenberg chain in a magnetic field as 
adapted from Ref. It models the magnetic excitations of 
the Hamiltonian Eq. ( p^ ) along q x = q y . The bar marks the 
width of the spectrum Aq for a given energy transfer. 

of the basal plane of the unit cell of Sr2Ru04, since the 
model is manifestly invariant under the symmetry trans- 
formations x «-> y and y — > —y, and since wave vectors 
are only defined modulus a reciprocal lattice vector we 
find gapless magnetic excitations with linear dispersion at 
Ql = (±[2s - 2k F ],±[% - 2fc F ])t » (±0.67r/o, ±0.67r/a) 
[compare Fig. |l| and Refs. 33 ]69|] . 

The result from conformed—field theory for any two 
point correlation functionE3o is valid for sufficiently 
small frequencies oj and momenta q = (q — |<?j|)<2 mea- 
sured with respect to the back scattering wave vector. 



Xm(?i U ) = 



v c gq 



2irT 



v c ft q 



2irT 



where 



Y(x/2-ik/2) 
T(l- x/2-ik/2) 



(33) 



(34) 



The value of the scaling dimension x and the excitation 
velocity v e s depend on the details of the system. The 
prefactor A x depends on the scaling dimension. Please 
refer to Ref. |7(] for details. 

The limits of the applicability of the result from con- 
formal field theory can be understood in the framework 
of studies of Heisenberg chains performed in Ref. [70]. It 
has been shown that the effective scaling dimension x is 
temperature dependent. At or above temperatures of the 
order of the excitation velocity, i.e., T > v e g we expect 
the scaling dimension to attain the non- interacting limit, 
i.e., x — ► 1. At energy transfers of the order of and above 
the excitation velocity lattice corrections become rele- 
vant. Similar arguments hold for the momentum trans- 
fer. At finite temperatures and finite energy transfer the 
effects combine and the r ange of validity of Eq. |33| can 
roughly be estimated as \/T 2 + u 2 + (v c $ q) 2 < v e g/2. 

The following experimental observations can be under- 
stood within the framework of the outlined analogies. 



100 



a} 



50 



O Experiment 
- Theory, jc=0.5, A =1.2 
Theory, x=0.65, A =1.1 




100 



200 



300 



T [K] 



FIG. 3. Plot of x m (0,6.2 meV) from Eq. (|3jJ) for 
x = 0.5 (broken line)r-and x — 0.65 (full line) in compari- 
son with experimentaEH results. Good agreement of the fits 
for fixed x can only be expected in the temperature interval 
2V m <T<.v cS . 



(i) The imaginary part of the magnetic correlation 
function at small energy transfer is strongly peaked 
at q i and flat elsewhere. This is in perfect qual- 
itative agreement with the magnetic structure factor 
S(q,w) ~ J«wr; m (q, w) determined via inelastic neutron 
scattering. 

(ii) The magnetic correlation function Eq. ( |33] ) is scale 
invariant. The scale invariance has been observed exper- 
imentally outside the Fermi liquid regimes and suggests 
values of 1/2 < x < 5/8 for Sr2Ru04. Note that these 
values of x describe a XXZ Heisenberg chain near the 
isotropic pointjld i.e., J > J z for in-plane (J) and out-of- 
plane (J 2 ) magnetic coupling ,L3 and are thus in quanti- 
tative agreement with the intermediate coupling regime 
assumed for Sr2Ru04 within this approach. 

(iii) Figure ^ shows the temperature dependence of 
lmx m (0, 6.2 meV) from Eq. ( j33| ) for x — 0.5 (bro- 
ken line) and-.:r = 0.65 (full line) in comparison with 
experimentaEa results. The fits are only valid outside 
the Fermi liquid regime, i.e., for T > 25 K. Together with 
the limitations of the applicability of Eq. ( p3[ ) discussed 
above the agreement of the fits for fixed x with the exper- 
imental data can only be expected in a small temperature 
interval of 2V m <T-^i v c g. Note the relatively large en- 
ergy transfer of u> — 6.2 meV = 73 K. Consequently a 
discrimination between x — 0.5 and x — 0.65 is not con- 
clusive. The prefactors A0.5 = 1.2 and A). 65 = 1.1 are 
of the correct order of magnitude. Neutron scattering 
results at lower energy transfer are desirable. 

(iv) The experimental results show a width of the 
magnetic r-aeaks which is only weakly temperature 
dependents as shown by the circles in Fig. [|. In the 
model presented here the finite width of the dynamic 
magnetic correlations follows out of the dispersion of the 
lower bounds of the excitation continuum as indicated by 
the bar near q ~ _ 2fc F in Fig. || with Aq ~ ^t^- 

The temperature dependence of the excitation veloc- 
ity can be estimated by considering the magnetic cor- 
relation function Im^ m (g, 6.2 meV) from Eq. (|33]) as a 
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FIG. 4. FuJi width at half maximum as determined by 
neutron scattering^! (circles) and as obtained from Eq. fel) 
with uj — 6.2 meV, x — 0.5, and v e g = 350 K (squares). The 
theoretical values for T > 75 K are less reliable since T and 
q become too large. The line is a guide to the eye. 



function of q and determining the full width at half max- 
imum (FWHM). The results are shown for x — 0.5 and 
w ff = 350 K as the squares in Fig. |J. Consistent with 
the expected limits of validity of Eq. ( |33| ) and the tem- 
perature dependence of the peak intensity discussed in 
paragraph (iii) the values of Aq are in good agreement 
with the experimental data only for T < 75 K. Note 
that for T > 75 K also the relevant values of q become 
large. For x = 0.65 the results only differ of the order 
of 10% and do not allow for any quantitative discrimi- 
nation. We conclude rj-tiat the magnetic energy scale is 
given by v c g ~ 10 2 K.O 

(v) The symbols in Fig. [s] show the energy dependence 
of lmx m (0,w) at T = 10.4 K from neutron scattering 
measurements. Open and full circles are from Ref. for 
energy and q scans, respectively. Squares from Ref. 
are scaled since no absolute scale is given. The data are in 
qualitative agreement with the presence of an excitation 
continuum. 

Since for T = 10.4 K the system is in the Fermi 
liquid regime the applicability of Eq. ( p3| ) is not ob- 
vious. The experimental data can be fitted with the 
(renormalized) non-interacting case, whercd x — 1 and 
lmx m (0, uj) ~ iQr tanh[w/(4r)] as shown by the full line 
in Fig. [5| The amplitude of the fit has been chosen as 
v c g = 120 K consistent with v c g ~ 10 2 K. The good 
agreement of the fit is likely to be accidental since for 
larger frequencies uj > 2 meV the system should gradu- 
ally cross over from Fermi liquid to conformal behavior. 
For lu (~ v c ff effects from the upper continuum limit be- 
come relevantd voiding the direct applicability of Eq. 
( |33| ) for lu > 10 meV. The broken line in Fig. || shows the 
result for the interacting case with x = 0.5 and Aq ^ = 1.2 
as determined under (iii) for comparison. 

Measurements of the low energy dynamical structure 
factor outside the Fermi liquid regime at 30 K < T < 50 
K are desirable to test the theory presented here in its 
range of applicability and allow conclusive comparison 
with perturbativeES results. 

(vi) The presence of quasi one-dimensional correlations 
along the system diagonals finds further experimental 
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5 10 
co [meV] 

FIG. 5. Plot of Imx m (0,o;)[T=io.4 k from Eq. @ for 
x = 1 (full line) and x — 0.65 (broken line) in comparisoo,with 
experimental^esults for frequency scans (open cirpfesEj and 
open square£3) and momentum scans (full circleaEB). Since 
at T ~ 10 K the system is in the Fermi liquid regime and 
for to > 2 me V only gradually crosses over to the conformally 
invariant regime the applicability of Eq. is not obvious 
(see text). 

support in the non-analytic .-angular dependence of the 
in-plane upper critical fields .Ej 

In conclusion the functional dependence of the dom- 
inant magnetic correlations anticipated from conformal 
field theory as given by Eq. (|||) describes the experi- 
mental data satisfactorily within the framework of the 
expected applicability of the theory. The consistency of 
the results suggests that we were able to extract a reliable 
energy scale for the effective magnetic correlations. 



B. Comparison with RPA 

The magnetic correlaticpSjJp— SrnR|iQ4 have been 
widely studied theoreticallyl^J'a'LErHcjO using pertur- 
bative approaches such as the random phase approxima- 
tion (RPA) . The perturbative approaches cannot account 
for the low dimensional quantum fluctuations and it is 
instructive to discuss the resulting limits of their appli- 
cability. To this end we have performed a RPA analysis 
of the magnetic structure factor. The interactio«-.is in- 
cluded in the dynamical correlation functions viaL-3 



1-17. 



V.\J 

(T.a' 



(35) 



The bare susceptibilities Xo' V (<Zi w ) are the Fourier trans- 
forms of the real time spin-spin correlation function 
—i6{S u q (t)S y _ q (0)) and are determined with respect to 



III. The 



the tight binding model discussed in Sec. 
spin operators S q (t) act on electrons with orbital index 
v,v' G {x,y,j}. In the absence of Hund's rule coupling 



the interactions are U^^, 



Uq and U 
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k z =0 




FIG. 6. Tr Im Xrpa(9' u ) convoluted with the experimen- 
tal q resolution. Parameters are u — 6.2 meV, Uo = 0.2 eV, 
Ui = 0.1 eV, [T'J, = t/£'J, = and dispersions Eqs. @ and 
(p^). Panel (a) k z — and (ib) fc z = 7r/c show tie sensitivity 
of the RPA approach to small changes in the parameters. 



For Ui ^ U2 the correlation function contributions to 
are also anisotropic in the spin Hilbert space. lB A 



tropic 
chO ii 



very recent approachO including Hund's rule coupling 
suggests the stabilization of a chiral magnetic state. The 
thermodynamic expectation values are determined via 
the interaction free Hamiltonian with tight binding bands 
as shown in Fig. [l| and given in Eqs. ([u]) and (|32|). 

Figure | shows Imx£p A (g,w) = TrImxRp A (9, ui) for 
u = 6.2 meV convoluted with the resolution 5q « 0.1%/a 
from neutron scattering experimentsEj (a) for k z = 
and (b) for k z — n/c. The interaction parameters here 
are U = 0.2 eV, U x =0.1 eV, and U^, = U^J, = 0. 
This neglects the hybridization of the d xy band with the 
d zx and c^—bands which has a quantitative effect on 
X Pi p A (<Z, but not so much on the total correlation 
function. The features discussed in Ref. [56] are repro- 
duced albeit with different weight. 

Panel (a) of Fig. |^ shows the correlations in the plane 
through the T point of the Brillouin zone (k z — 0) while 
the correlations in (b) lie in the_alane through the mid 
point of the line TZ (k z = 7r/c).Eil The difference of the 
two shows the sensitivity of the RPA approach to small 
changes in the parameters. 



The total correlation function in RPA in Fig. || clearly 
shows the structures of the bands dispersing in the 
plane. Many of these, structures have not been ob- 
served experimentally.E3 Moreover, the parameters have 
to be fine-limed in the RPA approach close to a phase 
transition.E3 Both effects can be understood as conse- 
quences of the underestimation of quantum fluctuations. 

Quantum fluctuations in low dimensions suppress long 
range order even if the value of the interaction is large 
enough, to give a finite temperature phase transition in 
RPA.O More specifically, for the one-dimensional case — 
which is discussed closely in the review by SolyomCa — 
the relevance of back and umklapp scattering terms 
is overestimated by perturbative approaches such as 
RPA. Instead, their relevance has to be determined non- 
perturbatively, namely via renormalization group (RG) 
studies. 

Consequently the two-dimensional RPA approach 
tends to underestimate the one-dimensional correlations 
since the renormalization of the excitation velocities can 
only be modeled indirectly by the interaction strength 
while the relative size of the two-dimensional features 
tends to be overestimated. Because of the large parame- 
ter space of the RPA approach and the sensitivity to de- 
tails in the band structure it is still possible to model the 
low temperature pagnetic structure factor in S^RuG^ 
rather accurately.E3 The weakness of the RPA approach 
becomes apparent through the fact that the description 
of different properties of ihe-.material requires different 
choices of parameter sets.E3 ! L3 

The present approach allows to include the RG results 
from the literatures since the back and umklapp scat- 
tering t erms are given explicitly by Eq. (|2l|) as discussed 
in Sec. [I A. The quantum fluctuations and interaction 



effects are included by the renormalization of parameters 
such as v e s and their impact becomes apparent thro ugh 



the small scaling dimension x discussed in Sec. [V A 



In conclusion the incommensurate magnetic fluctua- 
tions in Sr2Ru04 are best described via the quasi one- 
dimensional correlations-from Eq. (^3|) with a spectrum 
as sketched in Fig. |2|e3 The two-dimensional correla- 
tions yield an enhanced but relatively homogexLeous back- 
ground with only small additional structures £3 Note that 
the two-dimensional magnetic correlations beyond the di- 
agonals are given within the present approach by Eq. |l6| . 



C. Effective electronic masses 

In the following we will discuss effective electronic 
masses observed in the Fermi liquid regime of Sr2Ru04. 
Clearly, a system that shows scale invariant correlations 
down to 25 K must show a different renormalization of 
the Fermi liquid parameters than a system that is a Fermi 
liquid at all temperatures. Consequently it is quite nat- 
ural to estimate that difference phenomenologically by 
considering that there is a very strong reduction of the 
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effective excitation velocity of v-p /v c s ~ 20 — 60 in the 
magnetic sector along the diagonals of the basa l plane 
due to interaction effects as discussed in Sec. IV A[ In the 
Fermi liquid regime this reduction cannot be calculated 
direc tly via the results from conformal field theory (Sec. 
IV A ) because of the two-dimensional coupling V m (Sec. 
II A). It is yet reasonable to assume that the reduction 
of the excitation velocity in the magnetic channel along 
the diagonals remains much larger than the mere factor 
of 2 (see below) obtained in perturbativeEj approaches. 
With this assumption the different electronic masses can 
be modeled. 



1. Cyclotron mass 

The cyclotron mass m c determines the cyclotron fas- . 
quency and has been measured in dHvA experimentalxEa 
to be enhanced with respect to the bare (non-interacting) 
band mass mb on all three Fermi surfaces by the same 
amount, namely ra[ n, ''' 7 '/m[ ) a ^' 7 ' « 2. Since the cy- 
clotron motion does not involve magnetic excitations 
the cyclotron mass enhancement due to interactions 
is determined in the present approach by the quasi 
two-dimensional Hamiltonian of the charge channel Eq. 
( |j~5| ) in the subsystem of the d zx and d yz electrons. 
Perturbatiye. approaches with an appropriate choice of 



parametersEj also yield an electronic mass enhancement 



factor of 2 for all bands__ 
sistcnt with this result .LrE 



iRPES measurements are con- 



2. Specific heat 

The specific heat of an interacting system can be de- 
termined from the specific heat of the non interacting 
system via the renormaiization of the thermodynamic 
massc3 to* with respect to the bare band mass ox, 
equivalently, by the renormalized excitation velocitiesEJ 
Vp or v e s with respect to the bare band Fermi velocity 
v F . 

The specific heat of the d zx -d yz subsystem C z con- 
sists of the two-dimensional contributions from the spin 
and charge channels [Eqs. ( |l5| ) and (|l6|)] and the "one- 
dimensional" magnetic part along the diagonals [Sec. 



IV B|. The renormalized velocity of the two-dimensional 
contributions is given through the cyclotron mass en- 
hancement as v-p/vp 2D w 2. The quasi one-dimensional 
magnetic correlations with excitation velocity v c g are 
only present along the diagonals of the Brillouin zone. 
Consequently their contribution to the total specific heat 
must be weighed with respect to a system with one- 
dimensional magnetic correlations throughout the entire 
Brillouin zone. The normalized width of the magnetic 
peaks ^-Aq ~ 0.07 discussed in Sec. IV A allow for an 



estimate. Depending on whether one assumes quasi one- 
dimensional correlations in the vicinity of the positions 



of the incommensurate fluctuations q i only or along the 
whole diagonals weighing factors of 0.02 < w < 0.2 are 
reasonable. The weight of the two-dimensional magnetic 
contribution then is 1 — w. 

Together the contributions yield the specific heat of 
the d zx -d yz subsystem asc3 



a 



^7 = 7z 7o ~ y 



7o ■ (36) 



' J F,2D V F,2D v cS I 



The value of 7z « (m* /m ( h a) +2m* /m ( h fl) )/3 w 3.4 can be 
estimated from a weighed average of the thermodynamic 

masses m* /m]^ ss 3.1 and m*g/m^' » 3.5 measured via 
the dHvA effect Jl Then the "one-dimensional" magnetic 
contribution is determined as ip^E- k, 1.4 + w. Using the 
value of v c s — 350 K-jconsidered in Sec. IV A (iv) and the 
bare Fermi velocityEZl vp w 0.7 eV then gives w — 0.13. 

It must be pointed out that the two-dimensional cou- 
pling ~ V m <C Vm is likely to increase the effective ve- 
locity u e fi in the Fermi liquid regime through the reduc- 
tion of the correlation effects. This increase is compen- 
sated by an increase of the relevant phase space in the 
Brillouin zone determined by w. The phenomenological 
result that the strong low-dimensional magnetic correla- 
tions with v c g <C v F 2D enhance the specific heat beyond 
the value obtained via perturbative approachesEl remains 
unaltered. 

For non-interacting electrons in the two quasi one- 
dimensional bands under consideration the coefficient 
4.3 K ^ ol results in a specific heat contribution of 
that accounts for about 40% of the exper- 
imentally observed value of 7 to t = 40 ± 2 K ^ ol .E3~e3'e3 
Hydrostatic pressure increases the in-plane single par- 
ticle hopping which decreases the relative interaction 
strength and renders the d zx -d yz subsystem more two- 
dimensional. Consequently the "one-dimensional" mag- 
netic contribution in Eq. [36] is decreased yielding a 
natural explication for the observed reduction of the 
thermodynamic masses upon application of hydrostatic 
pressures within the present model. 

The thermodynamic mass of the d xy electrons is en- 
hanced by a factor pi m*/77i b 7 '' ~ 5.5 with respect to 
the bare band.jaaass.c3 The renormalized Ferpaa velocity 
from ARPESLl'L-3 or perturbative approaches^ only ac- 
counts for a factor ~ 2. The interaction between the 
d xy band and the d zx -d yz system accounts for a part of 
the missing enhancement through coupling to the one- 
dimensional correlations in the magnetic channel. An- 
other possible contribution comes from the proximity of 
the d xy band to the van Hove singularibijat the M point 
of the Brillouin zone.EJ Nesting effectsLliLZl yield an ad- 
ditional enhanced magnetic contribution to the specific 
heat. 



7o 

~T ~ 10 K%ol 
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3. Static susceptibility and Wilson ratio 

Following the argumentation of the specific heat the 
contributions to the uniform static magnetic susceptibil- 
ity of the d zx -d yz subsystem also consist of a two- and 
one-dimensional part. They are given with respect to 
the static magnetic susceptibility of the non-interacting 
system Xo a ^ 



i>f(1 — w) wv-p 



7 F,2D 



4-3 Xo • 



(37) 



All parameters have been fixed previously. Note that 
here the relative contribution of the "one-dimensional" 
subsystem is roughly twice that of the specific heat since 
7o includes both magnetic and charge degrees of freedom 
while Xo only accounts for the magnetic correlations. 

The enhancement of the susceptibility is in reasonable 
agreement with the relative spin- mass enhancement. p£ 
the a and (3 sheets measured via the dHvA effectEre 2 ] 
i a) « 3.7 and ml/m^ » 4.3. More- 



If A predicts the 



as m c 

over, since the model derived in Sec. 
quasi one-dimensional magnetic correlations only in the 
basal plane of the tetragonal lattice we expect a mag- 
netic mass enhancement that depends on the position k z 
on the FerrnLsurface with maxima at k Ztiaax {n) = 2irn/c 
as observedca experimentally. 

The ratio of the uniform static magnetic susceptibility 
and the specific heat coefficient has been determined ex- 
perimentally as Rw = (7r 2 x)/( 3 7tot) w lAEaaThe value 
of i?w,z ~ 1.3 obtained here for the d zx -d yz subsystem 
is in good agreement. The enhancement of the specific 
heat has indeed the same origin as the enhancement of 
the susceptibility as already concluded in Ref. [I] which 
can be identified in the d zx -d yz subsystem as the quasi 
one-dimensional magnetic correlations. 

An interesting experimental question that arises out 
of the discussion above is whether the photo electrons 
carry signatures of the enhanced magnetic correlations 
along the diagonals. To obtain an answer ARPES data 
for the a and (3 sheets need to be analyzed in detail along 
TX in compariso n wit h TM and MX. Using the value 
extracted in Sec. IV A (iv) the expected energy scale is 
VcB ~ 30 meV. 



V. CONCLUSIONS 



The dominantly one-dimensional kinetic energy of the 
electrons in the d zx and d yz orbitals allows to bosonize 
their Hamiltonian. In the presence of interaction this 
leads to an effective two-dimensional model. The degrees 
of freedom can be parameterized in terms of the four spin, 
charge, flavor, and spin-flavor fields. 

The presence of hybridization and corrections to the 
one-dimensional kinetic energy make the observation of 
properties of sliding Luttingcr liquids likely only at or 



above room temperature. This is consistent with the 
observed linear temperature dependence of the resistivity 
at high temperatures. 

In the magnetic sector described by the spin and spin- 
flavor fields the interaction at intermediate coupling leads 
to a quasi one-dimensional model along the diagonals of 
the basal plane of the Brillouin zone. The resulting spec- 
trum of elementary excitations accounts for the enhanced 
dynamical magnetic susceptibility at q i = (±2fcp,±2fcp) 
and the weak temperature dependence of the q width. 
The one-dimensional spectrum leads to a conformally in- 
variant formulation of the magnetic structure factor con- 
sistent with the observed scale invariance. The scaling 
dimension is consistent with the intermediate coupling 
regime. The observed excitation continuum and temper- 
ature dependence of the peak width are consistent with 
a magnetic energy scale of v c g ~ 10 2 K. The additional 
two-dimensional correlations are more homogeneous than 
predicted by RPA because of quantum fluctuations. 

The effective thermodynamic mass enhancement to- 
gether with the value of the specific heat coefficient is a 
superposition of two-dimensional effects as observed in 
ARPES and the "quasi one-dimensional" correlations of 
the magnetic channel of the d zx -d yz subsystem. The en- 
hancement of the static susceptibility has the same origin. 
The contribution of the quasi one-dimensional magnetic 
subsystem to the specific heat of the d zx -d yz subsystem 
is about 44% and about 70% to the static magnetic sus- 
ceptibility of the d zx -d yz subsystem. 
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